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c/3 ' The spatial variations of the order parameter, and the local density of states 

(LDOS) on the corner of s-wave or d x 2_ y 2-w&ve superconductors, as well as in 
superconductor-insulator-normal metal interfaces, are calculated self consis- 
tently using the Bogoliubov-deGennes formalism within the two dimensional 

CJ 

extended Hubbard model. The exact diagonalization method is used. Due 
■ to the suppression of the dominant d-wave order parameter, the extended s- 

wave order parameter is induced near the surface, that alternates its sign for 
the topmost sites at adjacent edges of the lattice and decays to zero in the 
bulk. The presence of surface roughness results into the appearance of the 
zero band conduction peak (ZBCP) near the corner surface which lacks from 



■ the predictions of the quasiclassical theory. 
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I. INTRODUCTION 



The determination of the order parameter symmetry has become one of the main as- 
pects in the research on high temperature superconductors MM. Tunneling conductance 
experiments report the existence of a zero band conduction peak (ZBCP) The ori- 

gin of the experimental ZBCP is explained in the context of zero energy states (ZES) || 
formed near the [110] surfaces of ci-wave superconductors These ZES do not appear 

for s-wave superconductors or near the [100] surface of ci-wave superconductors and are one 
of the features that characterize the ci-wave superconductors. 

The quasiclassical theory of superconductivity has been used to calculate the tunneling 
conductance in interfaces of unconventional superconductors with normal metals or ferro- 
magnets PHlO|- In the quasiclassical approximation the quasiparticles move in classical 
trajectories (external degrees of freedom) with internal degrees of freedom which are the 
spin and the particle-hole degrees of freedom. The orientation dependence of the spectra as 
well as the V line shape of the conductance curve are explained by the formation of bound 
states close to the interface due to the sign change of the pair potential that the transmitted 
quasiparticles experience. 

Moreover the concept of phase shift by tt of the order parameter in orthogonal directions 
in /c-space, which is equivalent to the sign change of the Josephson critical current, can 
be observed in corner junctions of anisotropic superconductors with conventional s-wave 
superconductors as a dip of the Fraunhofer pattern at zero magnetic field ||. It is an 
indication of ci-wave symmetry of the order parameter. The spontaneous flux modulation 
with surface orientation in such junctions has been calculated and can be used to distinguish 
the sub dominant components s or d xy that are induced at regions where the ci-wave order 
parameter is suppressed. 



The extended Hubbard model has been used to study single vortex structure [11,12 



time reversal symmetry breaking across twin boundaries |13| or near surfaces 0, the effect of 



disorder WM, and the effect of surface roughness H. In this paper the Bogoliubov-deGennes 
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equations are solved in a two dimensional square lattice within the context of an extended 
Hubbard model. The spatial variation of order parameter, and the local density of states, 
which in the limit of low transparency barrier converges to the tunneling conductance, are 
calculated for various types of surfaces and interfaces e.g. a corner surface, the interface of d- 
wave or s-wave superconductor along the [110] direction with normal metals. The evolution 
of the local density of states (LDOS) is studied as a function of the distance from the surface. 

It is seen that the extended s-wave order parameter is induced due to the suppression 
of the dominant rf-wave order parameter which alternates its sign for the topmost sites at 
adjacent edges of the lattice and decays to zero in the bulk. The LDOS is symmetric when 
(j, — and it becomes asymmetric when \i deviates from zero due to the breakdown of the 
electron-hole symmetry. We also investigate the effect of the surface roughness near the 
corner. In general surface roughness which in real samples is of atomic length scale modifies 
the properties of the quasiparticles since the coherence length of the cuprates is much smaller 
than the conventional s-wave superconductors. Our model treats the quasiparticle properties 
on atomic length scale and goes beyond the quasiclassical approximation. The presence of 
surface roughness results into the appearance of the ZBCP near the corner surface which 
lacks from the predictions of the quasiclassical theory. 

The rest of the paper is organized as following. In Sec. II we develop the model and 
discuss the formalism. In Sec. Ill we present the results for the corner of superconductor. 
In Sec. IV we present the s — I — n, d — I — n interfaces. In Sec. V, the effect of the surface 
roughness is considered. Finally, summary and discussions are presented in the last section. 

II. BDG EQUATIONS, WITHIN THE HUBBARD MODEL 

The Hamiltonian for the extended Hubbard model on a two dimensional square lattice 

is 

h = -t c lcja + n Yl n * + Pi 71 * 

<i,j>cr ia in 
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+ VqJ2 nil nil + y X n i° n j*' > 



(1) 



<ij><7(7 

where z, j are sites indices and the angle brackets indicate that the hopping is only to nearest 
neighbors, n ia = c\ a Ci a is the electron number operator in site i, fi is the chemical potential. 
Vo, V\ are on site and nearest-neighbor interaction strength. Negative values of V and V\ 
mean attractive interaction and positive values mean repulsive interaction. When V\ < 
the pairing interaction gives rise to the <i-wave superconductivity in a restricted parameter 



regime |T!|. To simulate the effect of the depletion of the carrier density at the surface or 
impurities the site-dependent impurity potential //(r^) is set to a sufficiently large value at 
the surface sites. This prohibit the electron tunneling over these sites. Within the mean 



field approximation Eq. ([[]) is reduced to the Bogoliubov deGennes equations (15 







u n {ri) 




' u n (n) 










I V n (ri) j 



(2) 



(3) 
(4) 



with 

iunin) = -tJ2 u n(ri + 5) + (//fa) + fi)u n (ri), 
5 

An„(rj) = A (rj)u n (ri) + ^ A s (ri)u n (ri + 5), 

6 

where the gap functions are defined by 

Aofc) = V < c^c^n) >, (5) 

A s (n) =V X < c T (r, + S) Cl (n) >, (6) 
where 5 = x, —x, y, —y. Equations (§) are subject to the self consistency requirements 

Ao(r,) = Vo E MnKiri) tanh {&\ , (7) 

V 



+ u n (ri)v*(ri + 5)) tanh 



(8) 



We start with the approximate initial conditions for the gap functions (0J1)- After exact 
diagonalization of Eq. (0) we obtain the u{rj) and v (r») and the eigenenergies e n . The 
quasiparticle amplitudes are then inserted into Eq. and new gap functions A (rj) and 
As(ri) are evaluated. We reinsert these quantities into Eq. (|,f§), and we proceed in the 
same way until we achieve selfconsistency, i.e., when the norm of the difference of A (rj) 
and As(ri) from their previous values is less than the desired accuracy. We then compute 



the d-wave and the extended s-wave gap functions given by the expressions fll 



A d { ri ) = ^[Az( n ) + A_ a (r;) - A § ( ri ) - A^(n)}, (9) 

AT{n) = \[A & (n) + A_*(r,) + A s (r,) + A_*(r,)]. (10) 
The number density at the ith site is given by 

n . = n ^ + Uil = J2 [\u n (ri)\ 2 f(e n ) + \v n (ri)\ 2 (l - /(e„))] , (11) 

n 

and the local density of states (LDOS) at the ith site is given by 

Pl (E) = -2^2 [\u n ( ri )\ 2 f'(E - e n ) + \v n ( ri )\ 2 f'(E + e n )] , (12) 

n 

where the factor 2 comes from the twofold spin degeneracy, / is the derivative of the Fermi 
function, 

= 1 it T\ I 1 ( 13 ) 

exp(e/k B T) + 1 



III. CORNER OF SUPERCONDUCTOR 

In this section we study the spatial variation of the order parameter close to the corner 
surface of a two dimensional square lattice seen in Fig. |1|. The different symmetries are 
introduced by varying the strength of the local and non-local pairing interaction constants 
Vq, Vl 
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We consider a two dimensional system of 30 x 30 sites, and we suppose fixed boundary 
conditions by setting the impurity potential fi 1 = lOOt at the surface. The temperature 
is ksT = O.lt. The constants Vq and V\ are 0.0 and — 2.5t respectively. This choice of 
parameters gives (i-wave superconductivity. The A^ order parameter is shown in Fig. |2|, 
for two different values of the chemical potential \l = and /i = t. For \x = the bulk value 
of the gap is A^ = 0.32t, while the transition temperature is calculated as /cbT c = 0.6t. It is 
much more suppressed for finite values of the chemical potential. As we can see the <i-wave 
order parameter is enhanced near the surface from its bulk value, and goes to zero at the 
surface atoms because the hopping to these sites is difficult due to the impurity barrier. In 
Fig. |3| we see the induced extended s-wave order parameter A e s xt , for two different values 
of the chemical potential \x = and \x = t. It is much more enhanced for greater values 
of the chemical potential. Near the surface A e s xt oscillates at an atomic scale and vanishes 
into the bulk region at a distance of few lattice sites. It reverses its sign on either side of 
the lattice edge and it is exactly zero in the diagonal direction. Next to the corner we see 
an enhancement from the edge value. It appears to have a <i-wave like structure just at 
the corner of the square lattice. This behavior is also seen near impurities [TJ|, and across 



twin boundaries [|T3|] using BdG equations within the extended Hubbard model in a two 
dimensional orthorhombic lattice. The explanation beyond this is the sign change of the 
d-wave order parameter across the [110] direction close to the corner. The quasiparticles 
feel the sign change of the pairing potential when they are reflected from the lattice edges 
near the corner. In Fig. |](a) we plot the local density of states for // = at points 
A, B, C, D seen in Fig. [I], along the diagonal of the lattice. The symmetry with respect 
to E = is due to the zero chemical potential. The LDOS is site depended and shows a 
complicated gap structure. Also no ZBCP has been observed in agreement with the results 
of the quasiclassical theory ||. As we move to the interior of the lattice the LDOS converges 
to the bulk density of states in a two dimensional square lattice. When [i deviates from zero 
as in Fig. |](b) the LDOS becomes asymmetric. However no ZBCP is formed. In Fig. |^(a) 
we plot the local density of states for /x = 0, at sites x = 1, 2, 3 seen in Fig. [[], along the 
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direction [100] perpendicular to the lattice edge. The LDOS has the VMike line shape, is 
symmetric with respect to E — due to the electron-hole symmetry, and has the minimum 
at E = which characterizes the d-wave symmetry of the order parameter. It is seen that 
the LDOS recovers its bulk value (long dashed line) as we move away from the surface. The 
symmetric form of the LDOS line shape is lost when the chemical potential deviates from 
zero, as seen in Fig. |(b) for jj, = t. 

To understand the effect of the different symmetry the local pairing interaction is set 
to the value Vq = —2.5t, while the non-local interaction is set to the value V\ = 0. This 
set of parameters gives rise to s-wave pairing. The spatial variation of the s-wave order 
parameter close to the corner of a two dimensional square lattice (see Fig. is seen in 
Fig. 0, for two different values of the chemical potential \x = and /x = t. For \x = as 
seen in Fig. |7](a) the order parameter increases monotonically (which is to be compared 
to the rf-wave case where this increase is non monotonous) to the bulk value at the sides 
of the lattice while at the corner an enhancement, relative to the bulk value is seen which 
disappears over a few sites to the interior of the lattice. For fi = t the order parameter 
increases non-monotonically even at the sides of the lattice as seen in Fig. 0(b). In Fig. ^a 
(b) we plot the local density of states at points A, B, C, D, for fi = 0(/i = t) seen in Fig. |5], 
along the diagonal direction of the lattice. The LDOS shows gap structure with [/-like line 
shape which characterizes the s-wave symmetry and is symmetric (asymmetric) when fi is 
(deviates from) zero . The LDOS is insensitive to the orientation of the surface and is site 
independent. Similar behavior is seen for the [100] surface seen in Fig. 

The absence of ZBCP in the LDOS can be explained from the results of the quasiclassical 
theory as follows. The condition for the formation of ZBCP is the change of sign of the 
quasiparticles during the scattering from the surface of the superconductor. In s-wave 
superconductors this sign change does not occur at surfaces or interfaces due to the isotropy 
of the pair potential. In anisotropic superconductors this sign change is possible for certain 
orientation of the surface. However for the corner surface, at the direction where the lobes 
of the o?-wave order parameter are at right angles to the surface, a typical trajectory of a 
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quasiparticle would consist of two subsequent reflections from the lattices edges, in non of 
which sign change of the order parameter occurs. Therefore the quasiparticle does not feel 
the sign change of the order parameter and no ZBCP is formed. For a corner where the lobes 
are not exactly at right angles to the surface the condition for the formation of Andreev 
bound states at the surface can occur and also the ZBCP. 

IV. S — I — N, D — I — N INTERFACES ALONG THE [110] DIRECTION 

We now discuss the effect of the interface on the order parameter and the local density of 
states for different symmetries. The interface is modeled by a line of sites along the diagonal 
of the lattice, y' direction, where the chemical potential is set to a value in accordance with 
the strength of the barrier we want to model. The value of the interaction strength in each 
part of the interface determines the particular system that we are considering. 

To understand the effect of the symmetry of the pair potential we consider first the 
conventional s-wave superconductor- insulator-normal metal (s — I — n) interface shown in 
Fig. |10|(a). The local interaction in region where x' < 0, (x' is the direction perpendicular 
to the interface) is Vq = —2.5t, and the strength of the barrier is /i 110 = lOOt. For \i = 
the superconducting gap is A s = 0.6t. We plot in Fig. |l0](b) the magnitude of the s-wave 
component A s of the superconducting order parameter as a function of x shown in Fig. |lO|(a) 
for two different values of the chemical potential, i.e., fi = and [i = t. It is suppressed near 
the interface and increases nonmonotonically to the bulk value at a few lattice sites. The 
enhancement at the topmost sites close to the interface is similar to the spatial variation 
of the ci-wave superconductor close to the [100] surface seen in Fig. 0. The bulk order 
parameter is suppressed when /i deviates from and also the spatial oscillations close to 
the interface are of reduced amplitude as seen in Fig. |10|(b). We also plot in Fig. |10|(c) the 
number density rii for the electrons for the two different values of the chemical potential, 
i.e., fi = and fi = t. For fi = the number density is unity (one electron per site) in the 
bulk and decays to zero at the interface. However for finite /i the number density is reduced. 
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In Fig. |TT|(b) we plot the local density of states (LDOS) at sites A, B, C, D, shown in 
Fig. [ll](a), for fi = 0. As expected the LDOS line shape is symmetric since /i = 0. The 
LDOS shows the gap structure and [/-like line shape. Also the LDOS is modified by the 
presence of the barrier, from it's bulk form. The LDOS is insensitive to the direction of the 
interface. For finite = 1 seen in Fig. |TT](c) the LDOS keeps its ?7-like line shape. However 
it becomes asymmetric due to the breakdown of the electron-hole symmetry. 

We consider then the unconventional <i-wave superconductor- insulator-normal metal (d — 
I — n) interface shown in Fig. 0(a) where V\ = — 2.5t, and /i 110 = lOOt. We plot in Fig. 
|l2|(b) the magnitude of the (i-wave component and the magnitude of the extended s-wave 
component A e s xt of the order parameter as a function of the direction x that is shown in Fig. 



12(a), for two different values of the chemical potential \i = and \i = t. It is seen that for 
/i = 0, A e s xt is not modified by the presence of the interface. In contrast A^ drops to zero 
at the interface, and increases monotonically into a few lattice sites to the bulk value. For 
fi = t the Ad is much more suppressed close to the interface while the induced A^-wave 
component is more enhanced. Due to symmetry arguments the relation A^(rj) = A^(r») 
holds for the pairing potentials and results into a monotonic variation in the <i-wave order 
parameter near the surface. This is different to the s — I—n case, and also to the [100] surface 
of the (i-wave superconductor seen in Fig. |2| where this increase is non-monotonous. In the 
latter case due to the translational invariance in the ^/-direction A_y(rj) = A^(rj) is satisfied. 
The spatial oscillations are due to the atomic scale oscillation of the A^r^), A_^(rj) that 
are completely neglected in the quasiclassical approximation. We also plot in Fig. [12|(c) the 
number density rii for the electrons for the two different values of the chemical potential, 
i.e., /i = and [i = t. It is similar to the case of the s — I — n- interface seen in Fig. |T0|(c) 

The main difference between the two symmetries appear in the LDOS as we can see 
in Fig. |13|(b) for the sites A,B,C,D, shown in Fig. 0(a), for \i = 0. As expected the 
LDOS line shape is symmetric since /z = 0. The ZBCP is formed and its height decreases 
exponentially as we move to the interior of the lattice along the direction perpendicular 
to the interface. However at site D no ZBCP is formed. The disappearance of ZBCP at 



D denotes that the zero energy states (ZES) wave functions have spatial variation close to 
the interface with nodes at specific sites. In this case site D corresponds to a node and 
therefore the ZBCP disappears. The ZBCP is explained in the context of zero energy states 
H formed near the [110] surfaces of d-wave superconductors due to the sign change that the 
quasiparticles experience in different directions in k space. However the absence of ZBCP for 
the [110] at specific sites is not predicted by the quasiclassical theory. For finite \x = 1 seen 
in Fig. |l3](c) the LDOS keeps its VMike line shape. However it becomes asymmetric due to 
the breakdown of the electron-hole symmetry. Also the ZBCP is reduced. We conclude that 
the order parameter as well as the local density of states are influenced by the orientation 



of the interface 116 



V. RESULT FOR THE SURFACE ROUGHNESS 

In the following we describe the effect of the surface roughness near the corner of the 
lattice. For the case of one step structure, for fi = 0, shown in Fig. 0(a) the LDOS 
shows ZBCP at points A,B but not in C as presented in Fig. 0(b). Moreover the ZBCP 
is suppressed compared to the case of flat [110] surface seen in Fig. [13]. The quasiclassical 
theory predicts that for that direction the ZBCP is maximum |§ since the one step structure 
corresponds to the a = 7r/4 (where a is the orientation of the surface). The suppression of 
the ZBCP and also the disappearance of ZBCP from specific sites is explained by the spatial 
variation of the ZES. It is seen that the wave functions of the ZES form standing waves that 
decay in the bulk. The sites A and B, that show ZBCP, correspond to an antinode while 
for the rest of the lattice sites the amplitude of the ZES is zero. For finite /i the ZBCP 
disappears from sites A and B due to the destructive interference of the ZES as shown in 
Fig. 0(c). In addition the overall line shape of the conductance curve is asymmetric due 



to the breaking of the electron-hole symmetry. In Fig. the spatial variation of the d- 
wave and extended s-wave order parameter is plotted for jj, — at sites close to the lattice 
corner. It is seen that the <i-wave order parameter is suppressed at the impurity site while 
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the extended s-wave order parameter is not much influenced. 

For the 1x2 step structure shown in Fig. 0(a) the LDOS, presented in Fig. |l6|(b) 
for fi = 0, at points A,B,C,D shows no ZBCP. The quasiclassical theory predicts a ZBCP 
since this geometry corresponds to a surface tilted from a = 0ora = 7r/2. The absence 
of ZBCP is explained by the destructive interference of the standing waves and also by the 
asymmetricity of the structure. For \x = t shown in Fig. [161(c) some tiny conductance peak 
recovers. 

The ZBCP at the topmost sites recovers for the geometry shown in Fig. |l7|(a) for sites 
A,C, for jj, = as presented in Fig. [13(b). However unlike the flat [110] surface no ZBCP is 
formed in site B due to the destructive interference of the standing waves. Moreover when 
finite chemical potential is introduced as shown in Fig. [13(c) the spatial distribution of the 
ZES is disturbed and for example the ZBCP appears even for sites, e.g., D, where normally 
for n = is absent. 



VI. CONCLUSIONS 

We calculated the LDOS and the order parameter of a two dimensional lattice of ci-wave 
superconductor within the extended Hubbard model, self consistently. The dominant order 
parameter decays monotonically for [110] interface and ci-wave order parameter symmetry, 
while non-monotonically for [100] interface, ci-wave or s-wave order parameter symmetry. 
The induced extended s-wave that decays to zero in the bulk changes sign at the topmost 
sites at either side of the lattice similarly to the case near impurities, and twin boundaries. 
The LDOS is symmetric when fi = and it becomes asymmetric when /i deviates from 
zero. The presence of surface roughness at the corner strongly modifies the quasiparticle 
properties near the corner. The ZBCP, which is absent for perfect corner consistently to the 
quasiclassical theory, appears when the roughness at atomic size is introduced due to the 
oscillatory form of the bound states that decay in the bulk. The last result is beyond the 
quasiclassical approximation. The sensitivity of the properties on the atomic scale has to 
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be taken into account for the correct interpretation of the experiments on corner junctions 
and in models that calculate the Josephson current phase relation, or the critical currents 
in those junctions taking into account the Andreev bound states. 



12 



REFERENCES 

[1] D.J. Scalapino, Phys. Rep. 250, 329 (1995). 

[2] D.J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 

[3] M. Covington, R. Scheuerer, K. Bloom, and L. Green, Appl. Phys. Lett. 68, 1717 (1996). 



[4] J. Lesueur, X. Grison, M. Aprili, and T. Kontos, |cond-mat/99092i~2 . 
[5] C.R. Hu, Phys. Rev. Lett, 72, 1526 (1994). 

[6] Y. Tanuma, Y. Tanaka, M. Yamashiro, and S. Kashiwaya, Phys. Rev . B 57, 7997 
(1998). 

[7] J.-X. Zhu, B. Friedman, C.S. Ting, Phys. Rev. B 59, 3353 (1999). 

[8] N. Stefanakis, J. Phys.: Condens. Matter 13, 1265 (2001). 

[9] N. Stefanakis, Phys. Rev. B to appear (2001). 
[10] N. Stefanakis, J. Phys.: Condens. Matter 13, 3643 (2001). 
[11] P.I. Soininen, C. Kallin, and A. J. Berlinsky, Phys. Rev. B 50, 13 883 (1994). 
[12] Y. Wang and A.H. MacDonald, Phys. Rev. B 52, 3876 (1995). 

[13] D.L. Feder, A. Beardsall, A. J. Berlinsky, and C. Kallin, Phys. Rev. B. 56, 5751 (1997). 
[14] M. Franz, C. Kallin, and A. J. Berlinsky, Phys. Rev. B 54, 6897 (1996). 
[15] P.G. de Gennes, Superconductivity of Metals and Alloys (Benjamin, New York, 1966). 
[16] A.M. Martin and J.F. Annett, |cond-mat/9811389 . 



13 



FIG. 1. The square lattice close to the corner. The chemical potential is set to fi 1 = lOOt at 
the shaded surface sites. 




FIG. 2. Spatial dependence of the d-wave order parameter close to the corner of a two dimen- 
sional square lattice, (a) fj, = 0, (b) /i = t. 
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FIG. 3. Spatial dependence of the extended s-wave Af* order parameter close to the corner of 
a two dimensional square lattice, (a) ji = 0, (b) fx = t. 
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FIG. 5. The local density of states at sites x = 1,2,3 shown in Fig. 1 along the [100] surface 
of the two dimensional square lattice, and the bulk density of states, (a) /i = 0, (b) /i = t. 



17 



FIG. 6. The square lattice close to the corner. The chemical potential is set to [i 1 = lOOt at 
the shaded surface sites. 
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FIG. 7. Spatial dependence of the s-wave order parameter close to the corner of a two dimen- 
sional square lattice, (a) \i = 0, (b) fi = t. 
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FIG. 9. The local density of states at sites x = 1,2,3 shown in Fig. 6 along the [100] surface 
of the two dimensional square lattice, and the bulk density of states. The pairing symmetry is s. 
(a) n = 0, (b) \i = t. 
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FIG. 10. (a) Spatial distribution of impurities, indicated as solid circles, corresponding to a 

s — I — n interface along the [110] surface. Also the labeling of the sites along the x direction is 

shown, (b) The magnitude of the s-wave component A s of the superconducting order parameter 
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as a function of x, for a s — I — n interface for fi = 0, t. The order parameter is calculated along 



the thick dashed line in direction x shown in (a), (c) The number density rii as a function of x 
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FIG. 11. (a) Spatial distribution of impurities, indicated as solid circles, corresponding to a 
s — I — n interface along the [110] surface, (b) The local density of states (LDOS) at A,B,C,D, 
shown in (a) of a s — I — n interface for fi = 0. (c) The local density of states (LDOS) at A, B, C, D, 
shown in (a) of a s — I — n interface for fi = t. 
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FIG. 12. (a) Spatial distribution of impurities, indicated as solid circles, corresponding to a 

d — I — n interface along the [110] surface. Also the labeling of the sites along the x direction 

is shown, (b) The magnitude of the d-wave component and the extended s-wave component 
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(Ag Xt ) of the superconducting order parameter as a function of x, for a d — I — n interface for 
u = 0. t. The order parameter is calculated alone the thick dashed line in direction x shown in (a). 
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FIG. 13. (a) Spatial distribution of impurities, indicated as solid circles, corresponding to a 
d — I — n interface along the [110] surface, (b) The local density of states (LDOS) at A,B,C,D, 
shown in (a) of a d — I — n interface for fi = 0. (c) The local density of states (LDOS) at A, B, C, D, 
shown in (a) of a d — I — n interface for fi = t. 
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FIG. 14. (a) Spatial distribution of impurities corresponding to the corner surface, with a 
lxl step structure. The chemical potential is set to \r = lOOt at the shaded surface sites, (b) 
The local density of states at the specified sites A, B, C, D, E shown in (a) for /i = 0. (c) The local 
density of states at the specified sites A, B, C, D, E shown in (a) for [x = t. 
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FIG. 15. Spatial dependence of the (a) <i-wave and (b) the extended s-wave order parameter 
for the specified geometry of Fig. 14(a). 
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FIG. 16. (a) Spatial distribution of impurities corresponding to the corner surface, with a 
1x2 step structure. The chemical potential is set to [i 1 = lOOt at the shaded surface sites, (b) 
The local density of states at the specified sites A,B,C,D shown in (a) for fi = 0. (c) The local 
density of states at the specified sites A, B, C, D shown in (a) for fj, = t. 
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FIG. 17. (a) Spatial distribution of impurities corresponding to the corner surface, with a 
[110] structure. The chemical potential is set to /j, 1 = lOOi at the shaded surface sites, (b) The 
local density of states at the specified sites A, B, C, D shown in (a) for fj, = 0. (c) The local density 
of states at the specified sites A, B, C, D shown in (a) for \x = t. 



28 



